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Circadian rhythms are known to be important drivers of human activity and the recent availabil¬ 
ity of electronic records of human behaviour has provided fine-grained data of temporal patterns 
of activity on a large scale. Further, questionnaire studies have identified important individual 
differences in circadian rhythms, with people broadly categorised into morning-like or evening-like 
individuals. However, little is known about the social aspects of these circadian rhythms, or how 
they vary across individuals. In this study we use a unique 18-month dataset that combines mo¬ 
bile phone calls and questionnaire data to examine individual differences in the daily rhythms of 
mobile phone activity. We demonstrate clear individual differences in daily patterns of phone calls, 
and show that these individual differences are persistent despite a high degree of turnover in the 
individuals’ social networks. Further, women’s calls were longer than men’s calls, especially during 
the evening and at night, and these calls were typically focused on a small number of emotionally 
intense relationships. These results demonstrate that individual differences in circadian rhythms 
are not just related to broad patterns of morningness and eveningness, but have a strong social 
component, in directing phone calls to specific individuals at specific times of day. 


I. INTRODUCTION 

Human activity follows a circadian rhythm that is re¬ 
flected at the psychological, physiological and biochemi¬ 
cal levels m- This rhythm is driven by endogenous cel¬ 
lular mechanisms, but may be modulated by exogenous 
factors. Circadian rhythms are in general synchronized to 
the day-night cycle. Within this cycle, there are differ¬ 
ences between individuals. Notably, there are morning 
and evening types, those who wake up early and those 
who prefer to sleep late. This may result from intrinsic 
differences in the circadian pacemaker circuit; morning¬ 
ness and eveningness have also associated with gender 
and personality traits mi. Circadian patterns can be 
also detected in statistics aggregated over large numbers 
of individuals. They have been observed in a wide range 
of phenomena: in human suicidal acts, times of exhibiting 
unethical behaviour, times of sexual activity, and times 
of heart attacks EES]. 

Aggregate-level daily rhythms are known to appear in 
electronic records of human activity, from mobility m 
to Wikipedia and OpenStreetMap edits mum], activity 
on Twitter |14l . and the number of mobile phone calls 
per hour [22 UM- In addition to the day-night cycle, 
these patterns are modulated by a number of endogenous 
factors such as the daily work schedule, commuting pat¬ 
terns, and the activity patterns of one’s social circles. In 
studies of daily rhythms inferred from electronic records, 
the focus has typically remained at the aggregate level 
(e.g. the rhythms of cities in [16]). 

Thus previous work in this area has either used ques¬ 


tionnaires to examine broad patterns of morningness and 
eveningness 0HS3, or used electronic records to exam¬ 
ine circadian rhythms in more detail, aggregated over 
many individuals. In this paper, we use a unique lon¬ 
gitudinal dataset pil ng that combines questionnaire 
data with mobile phone records, allowing us to take an 
individual-centric point of view, and study in detail the 
daily rhythms of mobile telephone calls of individuals. 
Our aim is to go beyond observing circadian rhythms 
in call frequency alone. To this end, we look at three 
questions: First, are there clear individual differences in 
daily call frequency patterns, and if so, how persistent 
are these given network turnover i.e. as network mem¬ 
bers change, does the daily call pattern also change? Sec¬ 
ond, are there daily patterns in other observables than 
call rates - specifically, the social aspects of call diversity 
and who is being called? Are calls made more randomly 
at certain times of day and directed at certain contacts 
at other times? Third, do the properties of callers and 
callees, such as gender, explain some features of the ob¬ 
served patterns? 

To this end, we employ a data set that comprises the 
exact times and recipients of all outgoing mobile phone 
calls of 24 individuals (’’egos” in the following) for 18 
months. This data set consists of 74,124 phone calls al¬ 
together. This set was originally collected for the purpose 
of studying social network turnover over time m dur¬ 
ing the study, the participants finish high school and go 
to work or university (often in another city), which gives 
rise to major turnover in their personal networks. The 
mobile phone call data set is accompanied by 3 surveys on 
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the contacts who were called (’’alters” in the following), 
including their gender and information on kinship. 

Together, these data allow us to study the daily 
rhythms of calls in terms of numbers and in terms of 
recipients (who is called and when). We find that in 
terms of call frequency at each hour of day, each indi¬ 
vidual has their distinct, persistent pattern. These daily 
patterns persist for individuals despite a high degree of 
network turnover, and thus appear to be characteristics 
of individual egos, rather than dependent on the identify 
of specific alters. Within these patterns, there are clear 
variations in the entropy of called alters, indicating that 
certain times of day (evening and night, typically) are re¬ 
served for calling specific alters, whereas at other times 
the recipients of calls are more diverse. For female egos, 
these variations are accompanied by strong variations in 
the average duration of calls that increases towards the 
night - these long calls are typically made to friends in¬ 
stead of family members. 


II. RESULTS 


ence between daily call patterns. For each ego, we cal¬ 
culate two different distances: self (d se it) and reference 
(d re f). The self-distance d se if for an individual i is the 
average JSD between the call patterns in (Ii,/ 2 ) and 
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tance measures the divergence of patterns of different 
egos in one time interval. For each time interval we 
calculate JSD between daily patterns of eras i and j: 
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on average lower values than d re f, meaning 
more similarity between an ego’s consecutive daily pat¬ 
terns than between the patterns of different egos in one 
interval. On average, for each ego, 87% ± 12% of ref¬ 
erence distances are higher than self-distances. Com¬ 
paring average values of distances over all egos we get 
(d re f) = 0.083 ± 0.28 while (d se if) = 0.05 ± 0.22 ((d se if) < 
(d re f) with t = 6.98 and p <C 1CD 6 , two-sample unequal 
variance t-test). To validate these results with another 
method, we have used the £ 2 norm (see Methods), and the 
results qualitatively agree with JSD: (d se if) = 0.11 ±0.02 
and (d re f) = 0.14 ± 0.03 with t = 6.11, p < 10 -5 . 


A. Persistence of individual daily call patterns 

We begin by computing the daily call patterns for all 24 
egos. The data time span is divided to three consecutive 
6 -month intervals 1 1 , I 2 , and I 3 . For each ego and each 6 - 
month interval, we compute the average fraction of calls 
placed at each hour of the day. Considering 6 -month in¬ 
tervals separately allows investigation the persistence of 
any observed differences: were specific features of indi¬ 
vidual patterns due to random fluctuations alone, they 
would not persist over all intervals. 

The resulting daily call patterns for 8 representative 
egos (4 male, 4 female) for all intervals are displayed in 
Fig. g Two features clearly stand out: First, while the 
call patterns of all egos follow the day-night cycle and 
calls at night are infrequent, there are significant differ¬ 
ences between individuals. As an example, the ego whose 
pattern is displayed in panel a) makes more calls in the 
morning than others, whereas for the ego of panel g) there 
are frequent calls at late hours. Second, it appears that 
each individual’s specific patterns are rather similar in 
all 6 -month intervals. Both observations hold for all 24 
egos. This persistence is noteworthy, since it is known 
that at the same time, the social networks of these in¬ 
dividuals undergo major turnover [T 8 j. Because of this, 
the observed persistence points towards intrinsic driving 
forces behind the daily patterns, as these do not strongly 
depend on an ego’s personal network composition. 

The persistence of individual daily patterns is con¬ 
firmed with a more detailed analysis. Here we use the 
approach of Ref. iJS] to show that the daily call pat¬ 
terns of an individual in different time intervals are more 
similar than the patterns of different egos within one 
time interval. We use the Jensen-Shannon divergence 
(JSD) (see Appendix for details) to measure the differ- 


B. Alter-specificity in call patterns 

We next turn to the question of where the individual 
daily patterns come from, and study the extent of a social 
component -that is, alter-specificity- in the call patterns. 
One can conceive of two extreme cases: 1) The patterns 
are entirely endogenous and the rate of call activity at 
each hour of the day is intrinsic to the ego. In this case, 
the called alters are picked at random (however, with a 
weight proportional to the time-averaged fraction of calls 
to each alter). 2) The patterns are alter-specific, that is, 
calls to certain alters are placed at certain hours, and 
the daily pattern is a superposition of the alter-specific 
patterns. 

To assess the extent of alter-specificity, we divide the 
day into 6 -hour time spans (night: 0AM-6AM, morning: 
6 AM 12AM, afternoon: 12AM-6PM, evening: 6 PM- 
12PM), and for each alter and each time span, com¬ 
pute the relative call entropies H re \. First, call entropies 
F/orig are calculated from the original data for each 6 -hour 
span. To get the relative entropies H le \, these values are 
then divided by average entropies (H ie {) calculated for 
a reference model where the times of calls to all alters 
are shuffled on a weekly basis for each ego (see Meth¬ 
ods). If for a given 6 -hour span iF rel < 1, calls to certain 
alters are emphasized within that time span, whereas if 
H re \ ft: 1, there is no alter-specificity. 

The relative entropies H re 1 for the same 8 individuals 
as in Fig. |Tj are displayed in Fig. [3j together with aver¬ 
aged relative entropy for all 24 individuals over all three 
intervals. The average relative entropy is at its highest 
in the afternoon, with ( H le \ ) ss 1, indicating large diver¬ 
sity of called alters. (H re 1 ) has its lowest point at night, 
when the number of calls is also low (see Fig. |T]) . This 




3 






FIG. 1. The daily call patterns of 8 individuals (a-h). The red lines denote the average fraction of calls placed at the 
corresponding hour for each of the three intervals 1 1 , I 2 , and I 3 . The black line is the average call pattern of all 24 individuals 
over all intervals. Areas shaded green show where an individual’s fraction of calls exceeds the average, while areas shaded red 
show where it falls below the average. 


indicates that the few calls made at night are typically 
directed to specific alters. As with the call frequency pat¬ 
terns, Fig. |T] clearly points out that the entropy patterns 
of different egos are different (compare, e.g., panels d and 
e). Likewise, each ego’s patterns appear fairly persistent; 
however, there is more variation here, especially in the 
morning and at night when the call frequency is low and 
the entropy measures are as a result noisy. 

We next focus on the specific alters behind the low- 
entropy times of day. For this, we first count the total 
number of calls by each ego to each alter in each interval, 
and rank the alters of each ego according to this num¬ 
ber.Alter ranks based on number of calls are known to 
reflect both the level of emotional closeness between ego 



FIG. 2. Histogram of d se if and d re f calculated for each ego. 
This plot shows the results for all egos and all time intervals. 


and alter (as indexed on a standard psychological 1-10 
emotional closeness scale), and the frequency face-to-face 
contact between ego and alter [THJ. Then, for each 6-hour 
interval (morning, afternoon, evening, and night) we cal¬ 
culate the fraction of calls directed at the two top-ranked 
alters. These fractions are shown in Fig. [4j again for the 
same individuals as in Fig. |T] On average, it appears 
that the fraction of calls to two top-ranked alters of each 
ego increases towards late hours and is often the highest 
at night, when there is in general only a small number 
of calls and low relative entropy. The high fractions in¬ 
dicate that decrease of entropy towards night typically 
comes from calls mainly to top-ranked alters; note that 
again, there is individual variation and although the top- 
alter fractions are often similar across intervals, in some 
cases, interval I\ behaves differently. This interval corre¬ 
sponds to the participants finishing high school and the 
following summer holidays, so differences in call behavior 
can be expected. 

Because Figs. [3] and [4] point towards a correlation be¬ 
tween low entropy and calls to top-ranked alters, we 
next quantify this as follows: As the baseline levels and 
slopes of Fig. [4] have a lot of variation, we take each ego 
and their relative entropies and fractions of calls to top- 
ranked alters at each 6-hour interval. Then, we compute 
the Pearson correlation coefficient between entropies and 
top-alter fractions for all egos. Out of the resulting 24 
correlation coefficients, 15 were significant withp < 0.05, 
with two positive coefficients and 13 negative averaging 
at r ~ —0.72. Thus for more than half of the egos, low 
entropy is clearly associated with a high fraction of calls 
to top alters, while for almost all the rest, no conclusive 
results can be drawn. 

Since there are alter-specific communication patterns 
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FIG. 3. The relative entropies for the same 8 individuals as in Figl, calculated for 6-hour intervals (M: morning 6AM-12AM, 
A: afternoon 12AM-6PM, E: evening 6PM-0AM, N: night 0AM-6AM). (o): interval I\, (□): interval I 2 , (o): interval I 3 . The 
black line indicates the average relative entropy for all 24 individuals over all three intervals. 



time of day 


FIG. 4. The fractions of calls to the two top-ranked alters for the same 8 individuals as in Figl, calculated for the same 6-hour 
intervals as in Fig3 (M, A, E, N). (o): interval Ji, (□): interval I 2 , (o)\ interval I 3 . 


and the nature of communication depends on the time 
of day, we also look at call durations at different times. 
Here, we use data on ego and alter attributes from the 
conducted surveys. Previous studies have looked at gen¬ 
der differences in talkativeness as well as differences in us¬ 
age of phones(both for landlines and mobile phones) [T5k 
E], using data from different countries and age groups. 
Most of the recent studies of talkativeness suggest that 
men and women are similar mm- However in most 
studies which compare phone usage difference between 
men and women, women have been reported to have 
longer calls [23.; 24]. The differences in phone usage of 
males and females have been linked to their different so¬ 
cial roles EMU- Here, we add two more dimensions and 
look at call durations at different times of the day, as well 
as durations of calls to different types of social links (kin 
or friend/acquaintance). 


Fig.[5]shows that overall, the average durations of calls 
by females are longer than those of calls by males, and 
that the difference largely depends on the time of day 
such that it increases towards the evening and is highest 
at night. A closer look shows that this difference arises 
mostly from calls to friends. Male and female call dura¬ 
tions to kin are fairly similar and do not depend much on 
the time of day. When the gender of the called alters is 
analysed (Fig. [6|, it is seen that by far the longest calls 
are by female egos to male alters at night; again the dif¬ 
ferences are the smallest in the afternoon, i.e. when all 
egos are typically in a similar social setting (at school, 
work, or university). The finding agrees with previous 
studies which suggest that females have different bond¬ 
ing strategies and use phones for different purposes com¬ 
pared to men (26J [28]. Since nighttime calls are often 
targeted at top-ranking alters (who typically are emo- 
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FIG. 5. Average duration of calls made by males and females 
to their kin, friends, and all social contacts. 



time of the day 

FIG. 6. Comparison of the average duration of calls to so¬ 
cial contacts of the same and opposite gender, separately for 
females and males. 

tionally close jTS]), and the egos are in their late teens 
and are possibly experiencing emotionally intense rela¬ 
tionships with their romantic partners, it is likely that 
these long calls often relate to romantic relationships. 
Since we have a relatively small sample of individuals (24 
total), one might think that the high values for call dura¬ 
tions for females in afternoons and nights might only be 
caused by one or few females who make very long calls 
in those hours. To rule out this possibility, for each indi¬ 
vidual we compare call durations made in the morning or 
afternoon with duration of calls made in the afternoon or 
night, using two-sanrple unequal variance t-test. We see 
that for 9 out of 12 females p-values for this test are less 
than 1CU 6 , whereas only for 2 males out of 12 we have 
such small p-values. 


III. DISCUSSION 

In contrast to conventional studies on daily patterns 
and circadian rhythms in social networks that focus on 


aggregates of very large numbers of individuals, here we 
focused on a small but rich sample that combined ques¬ 
tionnaire and mobile phone data in order to be able to ex¬ 
plore in much greater detail features that characterise in¬ 
dividuals circadian rhythms. Our focus has been on three 
specific issues, namely (1) whether there are individual- 
specific patterns of calling that mirror previously demon¬ 
strated individual patterns in the way individuals al¬ 
locate their social capital to their alters; (2) whether 
the individual-specific patterns of calling are persistent 
in light of network turnover and (3) whether there are 
gender differences in calling patterns. We show that in¬ 
dividuals do indeed have different daily patterns of call 
activity. These patterns vary beyond simple morning¬ 
ness/eveningness, as measured in questionnaire studies, 
and appear to be characteristic of the individual, in much 
the same way as their characteristic way of distributing 
their social capital among their alters [18]. Thus these 
individual patterns are persistent, in that the pattern of 
distributing calls across the day is consistent across the 
three time periods, despite the high degree of network 
turnover in the 18 months of the study, associated with 
leaving school and entering work or University mm- 
We also showed that there are striking gender differences 
in call duration pattern across the day: while women’s 
calls are generally longer than men’s calls, this was es¬ 
pecially true during the evening and at night. Evening 
calls to males and to friends by female egos were espe¬ 
cially long, and often involved calls to specific individuals, 
usually the top-ranked alters, who may be boyfriends. 

Given that humans naturally spend the night asleep, 
the tendency for calls to exhibit a striking diurnal peri¬ 
odicity is not, of itself, especially surprising, of course. 
However, in our sample, the vast majority of calls were 
made between midday and late evening, with the bulk of 
these occurring in the 6-hour slot between 12AM-6PM 
Fig- 0 It is notable that, despite our essentially diur¬ 
nal nature, rather few calls were made before midday. 
Most calls made by our cohort of subjects are likely to 
have been purely social rather than functional (i.e. work 
or leisure-activity related). Weissner [25] reported that 
certain types of conversations (notably story-telling and 
social conversations) are much more common during the 
evening than during the day among !Kung San hunter 
gatherers, with conversations involving economic matters 
or social criticism taking place mainly during daylight 
hours. In her sample, 81% of evening fireside conver¬ 
sations involved storytelling (relaying of adventures or 
experiences, especially in far off places, or tales about 
myths, social conventions and rituals, experiences dur¬ 
ing trance states or real life travels). Casual observation 
suggests that we find social events in the evening more 
engaging than ones held during the day, and that sto¬ 
ries (and especially stories of the supernatural) told at 
night have an especial frisson. Our data suggest that 
this nocturnal intensity spills over into conversations be¬ 
tween close friends (but, interestingly, much less often 
between close kin). Why this should be remains unclear, 
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and would clearly merit more detailed study. 

Within this broad pattern, the individual differences 
in the distribution of calling, and particularly the persis¬ 
tence of these individual differences in the light of social 
network turnover, are strongly suggestive of some kind 
of personality characteristic. It is possible that these dif¬ 
ferences in personal style simply reflect individual differ¬ 
ences in circadian pattern BH5I, and are a consequence 
of the fact that some individuals are more active in the 
morning and others more active in the evening. It is 
perhaps less likely that the calling patterns are due to 
differences in individual sleep/wake cycles, since the de¬ 
mands of the working day are likely to have required 
everyone to be active in the morning and even early ris¬ 
ers are unlikely to have gone to bed by 6PM (note that 
circumstances such as unemployment do have effects on 
daily cycles of individuals, see [30]). However, it could be 
that, physiologically, morning people are more likely to 
feel motivated to be socially engaged in the daytime and 
evening people more likely to be so in the evening. The 
fact that some individuals find the evening hours par¬ 
ticularly attractive, while others prefer the day, remains 
intriguing in this context and obviously merits more de¬ 
tailed investigation. 

Notwithstanding the fact that some individuals are 
night-oriented and others day-oriented, it seems that 
many (though not all) egos prefer to call certain alters 
at night. These are typically the one or two individu¬ 
als (mainly males and friends) that have special status 
for the ego (Fig. [4]). We know from our detailed ques¬ 
tionnaire data that the individuals who egos call most 
often are those to whom they are emotionally closest, 
and those they have the most frequent face-to-face con¬ 
tact with [IS]. It seems that this is especially charac¬ 
teristic of female egos, and much less so of male egos. 
Unlike women, men do not call either their girlfriends or 
their same-sex best friends for long chats in the evenings 
(even though their girlfriends may call them). This strik¬ 
ing sex difference in whom actively makes the effort to 
call is reminiscent of the finding reported by Palchykov et 
al. [21], for a very large cellphone dataset, that younger 
women (in particular) are much more proactive in calling 
their primary male contact than are men. As such, our 
finding speaks to the importance of female choice in hu¬ 
man mating strategies - that having made up their mind, 
women are typically much more focused in pursuing and 
investing in their relationships, and especially romantic 
relationships. This striking difference between the two 
sexes may reflect women’s more intensely social nature 
compared to men. The fact that these ’’special” calls are 
reserved for the evening reinforces the suggestion that 
the hours of darkness have a special quality for certain 
kinds of social interactions and social relationships. 

For both sexes, it was much less common to call kin 
during the evening. This would reinforce the claim that 
relationships with kin are less fragile than those with 
friends, and hence require less persistent and less spe¬ 
cial servicing m- Reserving calls to these individuals 


for times of the day when they are, or might seem to be, 
more intimate may reinforce the sense that the relation¬ 
ship is special. In effect, kin relationships come for free 
by virtue of the fact that they are kin and ego is embed¬ 
ded in a densely interconnected web of relationships with 
them, and therefore require less active maintenance. In 
contrast, the quality of friendships deteriorates rapidly 
(within months) in the absence of sufficiently frequent 
contact [d 1321123] ■ 

Our study was, by the standards of most online net¬ 
work analyses, very small scale, of course. However, its 
merit is that we have a complete record of the identity 
of all Alters called by our subjects, as well as their rela¬ 
tionships with those alters. We know that the individu¬ 
als they call most often are those with whom they have 
the emotionally closest relationships m, and those with 
whom they are most socially active P3- In this respect, 
we are able to add a layer of detail that is not usually 
possible in large scale network studies. 

A strength of our study was that it combined detailed 
mobile phone records with questionnaire data. Thus we 
have information on the nature of the relationship be¬ 
tween egos and the alters they are calling, in terms of gen¬ 
der, kinship and emotional intensity. Further, we know 
from previous findings that the number and duration of 
phone calls relates to the emotional intensity of the re¬ 
lationship, as well as the level of social activity pmii]. 
This dataset therefore allows for analysis of the social 
nature of circadian rhythms, rather than simply examin¬ 
ing aggregate analyses of mobile phone activity or broad 
scale questionnaire data. Due to the intensive, longitu¬ 
dinal data of the data collection, the sample size is rel¬ 
atively small. However, the nature of the data allows 
us to add a level of individual detail on identity of the 
callers and callees that is usually not possible in large 
scale social network studies 
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APPENDIX 

1. Our data and its use 

Our dataset includes 18 months of (outgoing) call and 
text records of the 24 individuals. In this study we have 
only used call records (both to mobile phones and land¬ 
lines). Participants filled out 3 questionnaires in the be¬ 
ginning (month 0), in the middle (month 9) and at the 
end (month 18), about people in their social network. 
They identified each contact (alter) as kin or non-kin, and 
provided all the different phone numbers that one contact 
might possibly have, and therefore these records are very 
comprehensive and do not miss a part of communica¬ 
tion because an alter has several numbers (with multiple 
phone providers) and/ or uses landline. They also pro¬ 
vided other information such as gender, how emotionally 
close they are to the person and the frequency of face-to- 
face contact with the person. The data on phone calls was 
obtained from the fully time-stamped, itemised monthly 
bills provided by arrangement with the provider, with 
the agreement of the subscriber. Participants received 
free mobile phone accounts for 18 months in return for 
taking part in the study. The participants, questionnaire 
and mobile phone datasets are described fully in previ¬ 
ous publications HU IH|. Anonymized, aggregated data 
is available online (see the SI of ehd- 


2. Calculating daily patterns 

To calculate daily patterns of each ego, we have taken 
data from all days in the time interval of interest and 
have allocated each call to a 6-hour time bin based on 
its time stamp. We then count total number of events 
of each hour and divide it by total number of events (of 


that ego) during the the time interval, to get the fraction 
of calls in that hour. In each time interval, we only used 
data of complete weeks in that time interval. 

3. Measuring similarity of patterns 

The Jensen-Shannon divergence is a measure of the 
difference of two probability distributions. It is a form of 
Kullback-Leibler divergence (KLD); unlike KLD, it works 
for probability distributions that contain zero-valued el¬ 
ements. The JSD for two probability discrete distri¬ 
butions Pi and P 2 follows the formula JSD{P\,P 2 ) = 
H(\Pi + \P 2 ) - WH^Pi) - H(P 2 )], where Pi = p z {a) 
and Pi(a) is the fraction of calls to each alter, and H is 
the Shannon entropy ( H(P) = — ]Tp(a) logp(a)). 

We have also used the / 2 -norm as a way to verify 
our results calculated using JSD. f 2 -norm is a similar¬ 
ity measure of two distributions, which is defined as: 
t 2 = VEPi(a)-P2(aW ■ 

4. Entropy patterns and relative entropies 

We calculate the call entropy for a given hour (or range 
of hours) as follows: first, the fraction of calls out of all 
calls to each alter a, p a , is counted for the specified hour 
(range of hours). Then, the call entropy for this hour 
(range of hours) is computed as H 0 r i g = — J/ Q p a logp a - 
In order to obtain the relative entropy, we repeatedly 
shuffle the original data as follows: for each week, the 
times and recipient alters of all calls are randomly shuf¬ 
fled. This reference model corresponds to a situation 
where the original call frequency pattern and the num¬ 
ber of calls to each alter are the same, but no prefer¬ 
ence is shown to any specific alter at any specific time. 
Then, for each shuffled set of data, we calculate call en¬ 
tropy similarly as for the original data, and average over 
N = 1,000 realizations to get the average reference en¬ 
tropy (H re f). Finally, the relative entropy is obtained 
as H re 1 = H OI i g /(H Te { }. The shuffling for the reference 
model is done on a weekly basis in order to minimize the 
effects of long-term dynamics, such as declining numbers 
of calls to alters, or alters appearing for the first time 
within the studied 6-month interval. 
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